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We present our recent lattice calculation with dynamical quarks using the overlap fermion formu- 
lation, which has exact chiral symmetry. It is possible to compare our data of meson mass and 
decay constant with the prediction from the chiral perturbation theory. From such comparison, 
we investigate the convergence property of the chiral expansion. For Nt = 2, we observe that the 
prediction to NLO does not converge at the scale of kaon mass. Based on this fact, we extend the 
analysis to the Nf = 2 + 1 case and carry out the extrapolation to the physical mass point using 
the NNLO formulae. 
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1. Introduction 



In lattice QCD, numerical simulations are carried out with quark masses given as simulation 
parameters. Since the quark mass around the physical value makes the cost of numerical simulation 
highly demanding, data are usually obtained at masses significantly heavier than those in nature. 
Results at the physical mass point are obtained by an extrapolation of the data points. It is therefore 
crucial for the accuracy of lattice calculation to make a reliable extrapolation in the function of 
quark masses or, equivalently, in the pseudo-scalar quark masses. The chiral perturbation theory 
(ChPT) gives a theoretical guide for this extrapolation 

ChPT is an effective theory constructed from QCD based on the chiral symmetry and its spon- 
taneous breaking. This theory describes the physics in the low-energy region p 2 m m- % where 
Nambu-Goldstone pions dominate the dynamics of the system. One of the characteristics of ChPT 

is that the Lagrangian is written in terms of an expansion in p 1 : Jzf = ££2 + ££\ H , where ££i„ 

contains interactions among mesons of momenta &{p 2n ). Based on ChPT, physical quantities are 
thus expanded in p 2 and nv\. For the quantities such as meson masses and decay constants, the 
leading order (LO) contribution from Jz?2 is corrected by the next-to-leading order (NLO) terms, 
which consist of one-loop effects with =5^ and tree-level insertions of 5£$. It is possible to add even 
higher order corrections from multi-loop level diagrams and higher order Lagrangeans. 

Often, lattice data are fitted using the NLO ChPT formula as it is the best known functional 
form for the quantity of interest. But, in many cases, chiral extrapolations are carried out in the 
mass region below the cut off A x = 2\[2iif n »1.2 GeV without checking the convergence of chiral 
expansion at NLO. For the kaon sector, in particular, the convergence at the kaon mass rriK ~ 500 
MeV is a non-trivial issue. One of the problems for the convergence test is that conventional 
fermion formalisms on the lattice explicitly break chiral symmetry or flavor symmetry then allow 
for inconsistency between numerical data and continuum theory. 

In this article, we present our study of the convergence property of ChPT using the numeri- 
cal simulation with dynamical overlap fermion [Q], with which chiral symmetry as well as flavor 
symmetry is exactly preserved on the lattice. We generated two sets of gauge configurations with 
different number of flavors Nf ■ = 2 ^ and 2 + 1 [Ql, on which a series of phenomenological quan- 
tities including the kaon B-parameter, Bk [@], have been calculated. We test the convergence of 
ChPT through the fit of light meson masses and decay constants, i.e. n?\ and f n for Nf = 2 and m- % , 
m 2 K , f n and for Nf = 2 + 1, and determine of the coupling constants of the ChPT Lagrangean, 
i.e. low-energy constants (LECs). We also obtain the fundamental physical quantities such as 
the chiral condensate, up and down quark masses and strange quark mass as a result of the chiral 
extrapolation. 

In the following section, we briefly explain how to obtain the data points on the lattices, which 
is common for the calculations with different number of flavors. To discuss the issue of conver- 
gence, we describe a test of ChPT performed for the Nf = 2 case [^] in Section ||. Based on this 
test, in Section [J we present the chiral extrapolation of the Nf ■ = 2 + 1 data by using the NNLO 
ChPT formulae. 
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2. Getting data points 

We refer [|], ||] for the details of the generation of the gauge configurations. For Nf = 2, we 
generate 10,000 HMC trajectories on a 16 3 x 32 lattice at six different sea quark masses m sea while, 
for Nf ■ = 2 + 1 , we generate 2,500 trajectories on a 1 6 3 x 48 lattice for ten combinations of up-down 
and strange sea quark masses, i.e. five m u /s times two m/s. 

For Nf = 2 (2+1), we calculate 50 (80) pairs of the lowest-lying eigenmodes on each gauge 
configuration and store them on the disks. These eigenmodes are used to construct the low-mode 
contribution to the quark propagators. The higher-mode contribution is obtained by conventional 
CG calculation with significantly smaller amount of machine time than the full CG calculation. 
Those eigenmodes are also used to replace the lower-mode contribution in the meson correlation 
functions by that averaged over the source location (low-mode averaging) [f7|, |p. We extract meson 
mass from the exponential decay of the time-separated correlation function of pseudo-scalar oper- 
ator (P(t)P(O)). The decay constant, which is defined by the matrix element of the axial-current 
operator srfn, is obtained simultaneously using the PC AC relation d^stf^ = 2m q P. 

Throughout the Monte Carlo updates for both N/s, the global topological charge of the gauge 
configurations is fixed to zero. This is necessary to avoid discontinuous change of the Dirac eigen- 
value, which is numerically too-expensive. The artifact due to fixing the topology is understood 
as a finite size effect [§] in addition to the conventional finite size effect. For the physical size of 
our lattice L w 1.7 fm, the finite size effect could be sizable. We calculate both kinds of finite size 
effect from the analytic formulae based on ChPT [ |Iol , [FT). In particular, for the effect of the fixed 
topology, we make use of the numerical data of the topological susceptibility determined on the 
same lattice configurations [|l2|] . 

In order to obtain the physical quark mass, we need to renormalize bare quark mass on the 
lattice as m^ ren ' = Z m mq . We obtain the renormalization factor Z m by calculating scalar and 
pseudo-scalar vertex functions in the momentum space in the Landau gauge and applying the 



RI/MOM scheme []13[]. In extracting Z m from the vertex functions, we control the contamina- 
tion from the spontaneous chiral symmetry breaking by using the the low-mode contribution to 
the chiral condensate JT4|]. Using the perturbative matching factor known to 4-loop level and the 
extrapolation to the chiral limit, i.e. m sea = for Nf = 2 and m uc i = m s . = for Nf = 2 + 1, we 
obtain the final results of Z m in the MS scheme at 2 GeV. 

In the rest of this article, it is understood that all data points are corrected by the finite size 
effects and quark masses are renormalized. 

3. Convergence of ChPT (Nf = 2) 

Determining the lattice scale by the Sommer scale with an input ro = 0.49 fm, we obtain oT x = 
1 .667 (17) GeV. It imply that our data points cover the pion mass region 290 MeV < m n < 750 MeV. 

In the framework of SU (2) ChPT which describes the data with Nf = 2, pion mass and decay 
constants are expanded in terms of x = ABm q / (4nf) 2 as 

m jr/ m <? = 2B{\ + 2*lnx) +c 3 x, (3.1) 
f% = f{l—x\nx) + c/\x (3.2) 
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Figure 2: Results of fit parameters B (top) and / (bottom) as functions of the upper limit of the fit range. In 
each panel, circle, square and diamonds are obtained with fit parameters x, x and t, . Results with ^ 2 /dof < 2 
are plotted. 



to NLO {i.e. one-loop level or &(xj), where B and / are the tree level LECs, and C3 and C4 are 
related to the one-loop level LECs Z3 and I4. At NLO, these expressions are unchanged when one 
replaces the expansion parameter x by x = 2m^/(4nf) 2 or £ = 2nv\j{A%f n ) 1 where rrr % and f K 
denote those at a finite quark mass. In other words, in a small enough pion mass region the three 
expansion parameters should describe the lattice data equally well. 

Three fit curves corresponding to x-fit, x-fit and <^ -fit for the three lightest pion mass points 
(m n < 450 MeV) are shown in Figure [I] as a function of m 2 K . For all fits, the horizontal axis is 
appropriately rescaled to give m 2 n using the obtained fit curves. From the plot we observe that the 
different expansion parameters describe the three lightest points equally well; the values of # 2 /dof 
are 0.30, 0.33 and 0.66 for x-, x- and ^ -fits, respectively. In each fit, the correlation between 
nv\jm q and f n for common sea quark mass is taken into account. Between the x- and x-fit, all 
of the resulting fit parameters are consistent. Among them, B and / are also consistent with the 
% -fit. This indicates that the NLO formulae successfully describes the data. In Figure ^[ results of 
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Figure 3: Comparison of the Nf = 2 results from the NLO fit and the NNLO fit with £, . Black pluses denote 
reference points from phenomenological estimations. 



B (upper panel) and / (bottom panel) from the different fits are plotted for various fitting range. 
As seen in the figure, the agreement among the different expansion prescriptions is lost when we 
extend the fit range to include the 4th lightest data point which corresponds to m K ~ 520 MeV. We, 
therefore, conclude that for these quantities the NLO ChPT may be safely applied only below as 
450 MeV. 

Another important observation from Figure [I] is that only the t, -fit reasonably describes the 
data beyond the fitted region. With the x- and x-fits the curvature due to the chiral logarithm is too 
strong to accommodate the heavier data points. In fact, values of the LECs with the x- and x-fits 
are more sensitive to the fit range than the -fit. This is because /, which is significantly smaller 
than f n of our data, enters in the definition of the expansion parameter. Qualitatively, by replacing 
m q and / by m„ and f„, higher loop effects in ChPT are effectively resummed and the convergence 
of the chiral expansion is improved. 

We then extend the analysis to include the NNLO terms [15]: 

ml/m q = 2B[l + ^ln^ + |(^ln^) 2 +(^-I(/~P h >' s + 16))^ln^+ C 3^(l-^ln^) + a^ 2 , 

(3.3) 

/* = /[l-§ln^ + |(^ln^) 2 + ^(/ phys + f )^ 2 ln^] +c 4 £ (1 - 5£ ln£) +/3 £ 2 . (3.4) 

Since we found that only the <^-fit reasonably describes the data beyond m n ~ 450 MeV, we perform 
the NNLO analysis using the t, -expansion. Although we input phenomenological estimate for the 
LEC P h y s , we find our fit result is insensitive to their uncertainties. We extract the LECs of ChPT, 
i.e. the decay constant in the chiral limit /, chiral condensate £ = Bf 2 /2, and the NLO LECs 
^hys _ _ C3 /g + \n(2\f2nf /m K +) 2 and /J hys = c\jf + \n(2^/2nf /m n +) 2 . For each quantity, a 
comparison of the results between the NLO and the NNLO fits is shown in Figure ||. In each panel, 
the results with 5 and 6 lightest data points are plotted for the NNLO fit. The correlated fits give 
# 2 /dof = 1.94 and 1.40, respectively. For the NLO fits, we plot results obtained with 4, 5 and 6 
points to show the stability of the fit. The # 2 /dof is less than 1.94. The results for these physical 
quantities are consistent within either the NLO or the NNLO fits. On the other hand, as seen for 
Z^ hys most prominently, there is a significant disagreement between NLO and NNLO. This is due 
to the large NNLO contributions to the terms which are proportional to C3 and C4, respectively. 
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Figure 4: Chiral extrapolation using the NNLO full SU (3) ChPT formulae for m\/m q (circles) and m\jm sli 
(squares) in the left panel and f n (circles) and fx (squares) in the right. 



We quote our results for the Nf = 2 calculation from the NNLO fit with all data points: / = 
111.7(3.5)(1.0)(+^)MeV,I M S(2 GeV) = [235.7(5.0)(2.0)(+ 1 ^) MeV] 3 , /f ys = 3.38(40) (24) (I 3 *), 
and /f ys = 4.12(35)(30)(± 3 (J), where m+ = 139.6 MeV. From the value at the neutral pion mass 
rrijfi = 135.0 MeV, we obtain the average up and down quark mass m uc i and the pion decay con- 
stant as mjf (2 GeV) = 4.452(81)(38)(+ 2 ° 7 ) MeV and f x = 119.6(3.0)(1.0)(+^) MeV. In these 
results, the first error is statistical, where the error of the renormalization constant is included in 
quadrature for E 1 / 3 and m uc i. The second error is systematic due to the truncation of the higher 
order corrections. For quantities carrying mass dimensions, the third error is from the ambiguity in 
the determination of r^. We estimate these errors from the difference of the results with our input 
ro = 0.49 fm and that with 0.465 fm. The third errors for if 1 ^ and /J hys reflect an ambiguity of 
choosing the renormalization scale of ChPT (4nf or 4nf K ). 

4. Results in the Nf = 2 + 1 simulation 

4.1 Fit to NNLO St/ (3) ChPT 

Since we found in the two-flavor calculation that the NNLO ChPT formulae can nicely fit our 
data even in the kaon mass region if one uses the E, -expansion, we apply the same strategy for our 
2 + 1-flavor analysis. As functions of £ n = 2m 2 K j{A%f n ) 2 and ^ = 2m 2 K / {An f n ) 2 , predictions from 
the SU (3) ChPT are expressed as 

m 2 Jm ud = 2B [l+M^^ K -,L r 4 ,L r 5 ,L r 6 ,L r 8 )] + af -g + ag + (4.1) 

m\lm sA = 2B [l+M K (^ K -,L r 4 ,L r 5 ,L r 6 ,L r 8 )]+af + of (4.2) 

h = fo[l+F^,^;Ll,L r 5 )] + P^C + ^-UK + Pi-^, (4.3) 

f K = fo [1 + F K fa, $k'AA)] + Pi ■ US* ~ k) + Pi ■ ~ (4.4) 

where m s d = \ {m s + m uc [) and af' and jS^' are NNLO unknown parameters. Functions M n , M K , 
F 71 and F K contain the NLO contributions from ^ and the NLO and NNLO loop contributions, 
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Figure 5: Chiral extrapolation using the NLO reduced SU (2) ChPT formulae. The organization is same as 
Figure 4. 



whose expressions are too involved to present here [16]. Among relevant SU(3) LECs, L^L^Lj 
and L' 7 which appear only in the NNLO contributions cannot be determined precisely. We intro- 
duce values L\ = (0.43 ±0.12) • 10~ 3 , L\ = (0.73 ±0.12) • 10~ 3 , L\ = (-2.53 ±0.37) • 10~ 3 and 



U-j = (-0.31 ± 0. 14) • 10~ 3 (defined at ju = 770 MeV) from a phenomenological estimate [|T7p and 



determine others L r 4 ,L r 5 ,L' 6 and Lg by a fit. Thus, the chiral extrapolation with Q4.1| )-(|4.4|) contains 
16 fit parameters in total. 

We fit m 2 n /m u d, m 2 K /m S( i, f% and fg simultaneously taking the correlation within the same sea 
quark mass (m uc i,m s ) into account. By using all data points, reasonable quality of the fit is obtained 
with ^ 2 /dof = 2.52. In this new study, we determine the lattice scale by the result of f n extrapolated 
to the physical point with the input f n = 130.0 MeV. As a result, we obtain a" 1 = 1.968(39) GeV 
and the pion mass covers the range of 340 MeV < m % < 870 MeV. Figure U shows all quantities in 
question as a function of m 2 . Different symbols correspond to the pion data (m 2 l /m lic j and f n ) and 
the kaon data (m^/ra^/ and fx) while the filled (open) symbols represent a fixed lighter (heavier) 
strange quark mass, which is accompanied by the solid (dashed) curves. 

Extrapolating the data to the physical point (^ phys ', ^i- phys '), which is determined with m % = 
135.0 MeV, m K = 495.0 MeV and f n = 130.0 MeV, we obtain preliminary results m^f (2 GeV) = 
3.64(12) MeV, m^(2 GeV) = 104.5(1.8) MeV, m s /m ud = 28.71(52), f K = 157.3(5.5) MeV and 
Ik/ fn = 1.210(12), where the errors are statistical only. 

In order to discuss the convergence property of ChPT as in the case of Nf = 2, we need to 
determine individual LECs with a high accuracy. However, with the data points obtained for two 
different strange quark masses, we have a limited constraint about the ^ dependence hence large 
errors for LECs. From the phenomenological side, it is advantageous to determine LECs along the 
line of this work because the results can be used as inputs in the calculation of different quantities 
including Bk and KIt, form factors. For these motivation, we are planning to extend the chiral 
extrapolation with more data points with m Ui / = m s . 
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Figure 6: Comparison of physical values between results from NNLO SU (3) analysis, NLO reduced SU (2) 
analysis and results obtained by UKQCD and RBC Collaborations [[I] 



4.2 Fit to the reduced 51/(2) ChPT to NLO 

As a check of the chiral extrapolation we carried out with the NNLO ChPT, we also study 
different fit ansatz. It is also possible to carry out the extrapolation to the physical point ^ phys ' by 
paying attention only to the dependence of the data on the up-down quark mass, or the pion mass. 
Integrating out the strange quark as a static heavy quark, one obtain an effective theory which 
respects a reduced SU (2) symmetry [18, 19, At NLO, the chiral expansion reads 



ml/m ud = 2B(l + i^ln^) +c 3 £„, 

m 2 K /m sd = 2B^+cP^, 

fx = /(l-£ t ln6 t )+C4&r, 

fK = f (K) (l-UMn)+cP^, 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



where we have LECs B^ K \ f^ K \ cf^ and in addition to the SU (2) LECs appeared in Section ||. 
In the present case, all LECs depend on strange quark mass. With the lightest three m ui i points, 
which are in the valid region of this framework, i.e. m ud <C m s for each fixed value of m s , we carry 
out the correlated fit for the quantities sharing the same mass point (m uc i,m s ). Figure || shows the 
fit curves obtained in this way. 



The fit results for each fixed m s are extrapolated to the physical strange quark mass m s 



(phys) 



which is determined by solving m 2 K /m s \ = (495.0MeV) 2 /m i . In FigureR we compare phys- 
ical results for m U( j, m s , fg and fn/fn fr° m the full NNLO SU (3) ChPT (circles), and from the 
NLO reduced SU (2) ChPT (squares from our analysis and diamonds from the similar analysis 
by RBC and UKQCD Collaborations [|^]). The agreement among different fitting prescription is 
encouraging. 



5. Summary 

We tested the convergence property of ChPT by comparing the analytic prediction with the 
lattice data obtained in the dynamical simulation with the overlap fermions. For Nf = 2, we carried 
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out the chiral fit to the NLO and NNLO formulae and compare three different expansion param- 
eters. We found that ChPT at NLO does not converge around the scale of kaon mass. It implies 
that one must take the NNLO effects into account to deal with the pion and kaon data points in an 
equal footing. In the Nf = 2 + 1 simulation, we fitted the data to the ChPT prediction to NNLO 
for the first time. The validity of the extrapolation to the physical mass point is checked with the 
results from the fit with the reduced SU (2) ChPT. We are planning to increase the data point with 
mud = m s to obtain the SU (3) LECs with high accuracy for a detailed study of the convergence 
property. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (Nos. 07-16 and 08-05 ). This work is supported in part by the Grant-in- Aid 
of the Ministry of Education (No. 20105005). 
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